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Abstract 

We study the imderljdng event in proton-antiproton collisions by examining the behavior of 
charged particles (transverse momentum px > 0.5 GeV/c, pseudorapidity |r|| < 1) produced in 
association with large transverse momentum jets (~2.2 fb ') or with Drell-Yan lepton-pairs {-2.1 
fb"') in the Z-boson mass region (70 < M(pair) < 110 GeV/c^) as measured by CDF at 1.96 TeV 
center-of-mass energy. We use the direction of the lepton-pair (in Drell-Yan production) or the 
leading jet (in high-px jet production) in each event to define three regions of ri-(|) space; toward, 
away, and transverse, where (j) is the azimuthal scattering angle. For Drell-Yan production 
(excluding the leptons) both the toward and transverse regions are very sensitive to the underlying 
event. In high-px jet production the transverse region is very sensitive to the underlying event and 
is separated into a MAX and MIN transverse region, which helps separate the hard component 
(initial and final-state radiation) from the beam-beam remnant and multiple parton interaction 
components of the scattering. The data are corrected to the particle level to remove detector 
effects and are then compared with several QCD Monte-Carlo models. The goal of this analysis is 
to provide data that can be used to test and improve the QCD Monte-Carlo models of the 
imderlying event that are used to simulate hadron-hadron collisions. 

I. INTRODUCTION 

In order to find physics beyond the Standard Model at a hadron-hadron collider, it is 
essential to have Monte-Carlo models that accurately simulate QCD hard-scattering events. To 
do this one must not only have a good model of the hard-scattering part of the process, but also 
of the beam-beam remnants (BBR) and the multiple parton interactions (MPI). The underlying 
event consists of the BBR plus MPI and is an unavoidable background to most collider 
observables. A good understanding of the underlying event will lead to more precise 
measurements at the Tevatron and the Large Hadron Collider (LHC). The goal of this analysis is 
to provide data that can be used to test and improve the QCD Monte Carlo models of the 
underljdng event. 

Figure 1 . 1 illustrates the way the QCD Monte-Carlo models simulate a proton-antiproton 
collision in which a hard 2-to-2 parton scattering with transverse momentum, pT(hard), has 
occurred. The resulting event contains particles that originate from the two outgoing partons 
(plus initial and final-state radiation) and particles that come from the breakup of the proton and 
antiproton. The beam-beam remnants are what is left over after a parton is knocked out of each 
of the initial two beam hadrons. They are one of the reasons why hadron-hadron collisions are 
more complicated than electron-positron annihilations. For the QCD Monte-Carlo models the 
beam-beam remnants are an important component of the underlying event. In addition to the 
hard 2-to-2 parton-parton scattering and the beam-beam remnants, sometimes there are 
additional semi-hard 2-to-2 parton-parton scatterings (MPI) that contribute particles to the 
underlying event. However, on an event-by-event basis these two components cannot be 
uniquely separated from particles that come from the initial and final-state radiation. Hence, a 
study of the underlying event inevitably involves a study of the BBR plus MPI plus initial and 
final-state radiation. 

As shown in Fig. 1.2, Drell-Yan lepton-pair production provides an excellent place to 
study the underlying event. Here one studies the outgoing charged particles (excluding the 
lepton pair) as a fimction of the lepton-pair invariant mass and as a fimction of the lepton-pair 
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transverse momentum. Unlike high-px jet production, for lepton-pair production there is no 
final-state gluon radiation. 




Fig. 1.1. Illustration of the way QCD Monte-Carlo models simulate a proton-antiproton collision in which a hard 2-to-2 parton 
scattering with transverse momentum, pT(hard), has occurred. The hard-scattering component of the event consists of particles 
that result from the hadronization of the two outgoing partons (i.e. the initial two jets) plus the particles that arise from initial and 
final state radiation (i.e. multi-jets). The underlying event consists of particles that arise from the beam-beam remnants and from 
multiple parton interactions. 




Fig. 1.2. Illustration of the way QCD Monte-Carlo models simulate Drell-Yan lepton-pair production. The hard-scattering 
component of the event consists of the two outgoing leptons plus particles that result from initial-state radiation. The underlying 
event consists of particles that arise from the beam-beam remnants and from multiple parton interactions. 

Hard-scattering collider jet events have a distinct topology. A typical hard-scattering event 
consists of a collection (or burst) of hadrons traveling roughly in the direction of the initial two 
beam particles and two collections of hadrons Gets) with large transverse momentum. The two 
large transverse momentum jets are roughly back-to-back in azimuthal angle, ((). One can use the 
topological structure of hadron-hadron collisions to study the underlying event. We use the 
direction of the leading (highest px) jet in each event to define four regions of ri-cj) space, where 
r\ is the pseudorapidity. The pseudorapidity r\ = -log(tan(6cm/2)), where 6cm is the center-of-mass 
polar scattering angle and (j) is the azimuthal angle of outgoing charged particles. As illustrated 
in Fig. 1.3, the direction of the leading jet, jet#l, in high-pxjet production or the direction of the 
lepton-pair in Drell-Yan production is used to define correlations in the azimuthal angle, A(j). 
The angle A(j) = (() - (()jet#i (A(j) = (j) - (j)pair) is the relative azimuthal angle between a charged 
particle and the direction of jet#l (lepton-pair). The toward region is defined by |A(j) | < 60° and 
|ri| < 1, while the away region is |A(j) | > 120° and |ri| < 1. The two transverse regions 60° < -A(() < 
120°, |ti| < 1 and 60° < A(() < 120°, |ti| < 1 are referred to as transverse 1 and transverse 2. The 
overall transverse region corresponds to combining the transverse- 1 and transverse-2 regions. In 
high-px jet production, the toward and away regions receive large contributions from the 
outgoing high-px jets, while the transverse region is perpendicular to the plane of the hard 2-to-2 
scattering and is therefore very sensitive to the underlying event. For Drell-Yan production both 
the toward and the transverse region are very sensitive to the underlying event, while the away 
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region receives large contributions from the away-side jet from the subprocesses: 

q + q ^ rr + g , q + g ^ l^T + q , q + g ^ l^T + q . 



Jet #1 Direction Lepton-Pair Direction 




Fig. 1.3. Illustration of coiTelations in azimuthal angle A(|) relative to (left) the direction of the leading jet (highest pj jet) in the 
event, jet#l, in high-pj jet production or (right) the direction of the lepton-pair in Drell-Yan production. The angle A(|) = (|) - 
(|)jet#l (A(j) = (|) - (|)pair) is the relative azimuthal angle between charged particles and the direction of jet#l (lepton-pair). The 
toward region is defined by |A(|) [ < 60° and [r|[ < 1, while the away region is |A(|) [ > 120° and |r|[ < 1. The two transverse regions 
60° < -A(|) < 120°, [r|[ < 1 and 60° < A(|) < 120°, [rjj < 1 are referred to as transverse 1 and transverse 2. Each of the two transverse 
regions have an area in r\-^ space of AriA(|) = 4n/6. The overall transverse region corresponds to combining the transverse- 1 and 
trans verse-2 regions. The transMAX (transMIN) refer to the transverse region (transverse- 1 or trans verse-2) containing the 
largest (smallest) number of charged particles or to the region containing the largest (smallest) scalar px sum of charged particles 

Table 1.1. Observables examined in this analysis as they are defined at the particle level and the 
detector level. Charged tracks are considered good if they pass the track selection criterion given in 
Section 111(5). The mean charged-particle <Pt> is constructed on an event -by-event basis and then 
averaged over the events. For the average px and the PTmax, we require that there is at least one 



charged particle present. Particles are considered stable if ex > 10 mm (K^, A, Z, S, and Q are kept 
stable) . 



Observable 


Particle Level 


Detector level 


dN/dr|d(|) 


Number of stable charged particles 
per unit rj-cj) 
(px > 0.5 GeV/c, |ri|< 1) 


Number of good tracks 

per unit rj-cj) 
(px > 0.5 GeV/c, |ri|< 1) 


dPT/drid(t) 


Scalar px sum of stable charged 
particles per unit r\-^ 
(px > 0.5 GeV/c, |ri|< 1) 


Scalar px sum of good tracks 
per unit r|-(j) 
(px>0.5 GeV/c, |ri|<l) 


<px> 


Average px of stable charged particles 
(px > 0.5 GeV/c, |ri|< 1) 
Require at least 1 charged particle 


Average px of good tracks 
(px > 0.5 GeV/c, |ri|< 1) 
Require at least 1 good track 


PTmax 


Maximum px stable charged particle 

(px > 0.5 GeV/c, |ri|< 1) 
Require at least 1 charged particle 


Maximum px good charged tracks 
(px > 0.5 GeV/c, |ri|< 1) 
Require at least 1 good track 


Jet 


MidPoint algorithm R = 0.7 f,nerge = 
0.75 apphed to stable particles 


MidPoint algorithm R = 0.7 fmeige ^ 
0.75 applied to calorimeter cells 



We study charged particles in the range pi > 0.5 GeV/c and |ri| < 1 in the toward, away 
and transverse regions. For high-px jet production, we require that the leading jet in the event be 
in the region |ri(jet#l)| < 2 (referred to as leading-jet events). The jets are constructed using the 
MidPoint algorithm (R = 0.7, fmerge = 0.75), where R is the jet radius and fmerge is the jet splitting 
and merging fraction [1]. For Drell-Yan production we require that the invariant mass of the 
lepton-pair be in the mass region of the Z-boson, 70 < M(pair) < 1 10 GeV/c , with |ri(pair)| < 6 
(referred to as Drell-Yan events). For both leading-jet and Drell-Yan events we define MAX and 
MIN transverse regions (transMAX and transMIN) [2], where MAX (MIN) refers to the 
transverse region (transverse- 1 or transverse-2) containing the largest (smallest) number of 
charged particles or to the region containing the largest (smallest) scalar px sum of charged 



Page 6 of 26 



particles. For events with large initial or final-state radiation the transMAX region will usually 
contain the third jet in high-px jet production or the second jet in Drell-Yan production while 
both the transMAX and transMIN regions receive contributions fi-om the beam-beam remnants. 
Thus, the transMIN region is very sensitive to the beam-beam remnants, while the event-by- 
event difference between transMAX and transMIN is very sensitive to initial and final-state 
radiation (transDIF = transMAX - transMIN). 

Table 1 . 1 shows the observables that are considered in this analysis as they are defined at 
the particle level and detector level. The detector level corresponds to the tracks passing good- 
track criteria and the particle level corresponds to true charged particles in the event. The 
particle level can be compared directly with the QCD Monte-Carlo models at the generator level. 
Since we will be studjdng regions in r\-^ space with different areas, we construct densities by 
dividing by the area. For example, the number density dN/dridcj), corresponds the number of 
charged particles per unit ri-(j) and the PTsum density dPT/drid(j), corresponds the charged-particle 
scalar-px sum per unit rj-c]). 

A discussion of the QCD Monte-Carlo Model is presented in Section II. In Section III we 
discuss the data selection, track cuts, and the method we use to correct the data to the particle 
level. Section IV contains the results for leading-jet and Drell-Yan events and comparisons with 
the QCD Monte-Carlo models. Section V is reserved for the summary and conclusions. 

II. QCD Monte-Carlo Models 

QCD Monte-Carlo generators such as PYTHIA [3] have parameters which may be adjusted to 
control the behavior of their event modeling. A specified set of these parameters that has been 
adjusted to better fit some aspects of the data is referred to as a tune. PYTHIA Tune A was 
determined by fitting the CDF Run 1 underlying event data [4]. Later it was noticed that Tune A 
does not fit the CDF Run 1 Z-boson px distribution very well [5]. PYTHIA Tune AW fits the Z- 
boson Px distribution as well as the underlying event at the Tevatron [6]. For leading-jet 
production Tune A and Tune AW are nearly identical. Table 2.1 shows the parameters for 
several tunes for PYTHIA version 6.2. PYTHIA Tune DW is very similar to Tune AW except the 
setting of one PYTHIA parameter PARP(67) = 2.5, which is the preferred value determined by the 
D0 Collaboration in fitting their dijet A(j) distribution [7]. PARP(67) sets the high-px scale for 
initial-state radiation in PYTHIA. It determines the maximal parton virtuality allowed in time-like 
showers. Tune DW and Tune DWT are identical at 1.96 TeV (the reference point), but Tune 
DW and DWT extrapolate differently to the LHC. Tune DWT uses the PYTHIA default value for 
energy dependence of the MPI cut-off (PARP(90) = 0.16), which is the value used in the ATLAS 
PYTHIA tune [8]. Tune DWT produces more activity in the underlying event at the LHC than 
does Tune DW, but predicts less activity than Tune DW in the underlying event at energies 
below 1.96 TeV. Tune DW uses the Tune A value of PARP(90) = 0.25, which was determined 
by comparing the Run 1 data at 1.8 TeV with the CDF underlying event measurements at 630 
GeV [9]. The amount of MPI and hence the tunings depend on the choice of the parton 
distribution functions. All these tunes use the CTEQ5L [10] parton distribution functions. 

The first 9 parameters in Table 2.1 tune the MPI. PARP(62), PARP(64), and PARP(67) tune 
the initial-state radiation and the last three parameters set the intrinsic transverse momentum of 
the partons within the incoming proton and antiproton. 
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Table 2.1. Parameters for several PYTHIA 6.2 tunes. Tune A is the CDF Run 1 underlying-event tune. Tune 
AW and DW are CDF Run 2 tunes whieh fit the existing Run 2 underlying event data and fit the Run 1 Z- 
boson Pt distribution. The ATLAS Tune is the tune used in the ATLAS TDR [8]. Tune DWT uses the 
ATLAS energy dependence for the MPI, PARP(90). The first 9 parameters tune the multiple parton 
interactions. PARP(62), PARP(64), and PARP(67) tune the initial-state radiation and the last three 
parameters set the intrinsic k-p of the partons within the incoming proton and antiproton. 



Parameter 


Description 


Tune 

A 


Tune 
AW 


Tune 
DW 


Tune 
DWT 


ATLAS 


PDF 


Parton Distribution 
Functions 


CTEQ5L 


CTEQ3L 


CTEQ3L 


CTEQ3L 


CTEQ3L 


MSTP(81) 


MPI On 


1 


1 


1 


1 


1 


MSTP(82) 


Double Gaussian 
Matter Distribution 


4 


4 


4 


4 


4 


PARP(82) 


MPI Cut-Off 


2.0 


2.0 


1.9 


1.9409 


1.8 


jr . V i\ 1 yoj J 


Fraction of matter 
within core 


0.5 


5 


5 


5 


5 


PARP(84) 


Core Radius 


0.4 


0.4 


0.4 


0.4 


0.5 


PARP(85) 


Color Connections 


0.9 


0.9 


1.0 


1.0 


0.33 


PARP(86) 


Color Connections 


0.95 


0.95 


1.0 


1.0 


0.66 




Reference Energy 


1800 


1800 


1800 


1960 


1000 


PARP(90) 


MPI Energy 
Dependence 


0.25 


0.25 


0.25 


0.16 


0.16 


PARP(62) 


Initial-state radiation 
Cut-Off 


1.0 


1.25 


1.25 


1.25 


1.0 


PARP(64) 


Soft Initial-State 
Radiation Scale 


1.0 


0.2 


0.2 


0.2 


1.0 


PARP(67) 


Hard Initial-State 
Radiation Scale 


4.0 


4.0 


2.5 


2.5 


1.0 


MSTP(91) 


Gaussian Intrinsic kj 


1 


1 


1 


1 


1 


PARP(91) 


Intrinsic Gaussian 
Width, o 


1.0 


2.1 


2.1 


2.1 


1.0 


PARP(93) 


Intrinsic kj Upper 
Cut-Off 


5.0 


15.0 


15.0 


15.0 


5.0 



Table 2.2. The computed value of the multiple parton scattering cross section for the various PYTHIA 6.2 
tunes. 



Tune 


o(MPI) 
at 1.96 TeV 


ct(MPI) 
at 14 TeV 


A, AW 


309.7 mb 


484.0 mb 


DW 


351.7 mb 


549.2 mb 


DWT 


351.7 mb 


829.1 mb 


ATLAS 


324.3 mb 


768.0 mb 



Table 2.2 shows the computed value of the multiple parton scattering cross section for the 
various tunes. The multiple parton scattering cross section (divided by the total inelastic cross 
section at the center-of-mass energies of 1.96 and 14 Tev, respectively) determines the average 
number of multiple parton collisions per event. The MPI cross section is the same for proton- 
proton and proton-antiproton collisions. 

HERWIG [11] is a QCD Monte-Carlo generator similar to PYTHIA except HERWIG employs a 
cluster fragmentation model while PYTHIA uses a string fragmentation approach. In addition, 
gluon radiation is modeled differently by the two generators. Also, HERWIG does not include 
MPI in the underlying event. In HERWIG the underlying event arises solely from the BBR. 
JIMMY [12] is a multiple parton interaction model which can be added to HERWIG to improve 
agreement with the underlying event observables. To compare with the Drell-Yan data we have 
constructed a HERWIG tune (with JIMMY MPI) with JMUEO = 1, PTJIM = 3.6 GeV/c, 
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JMRAD(73) = 1.8, and JMRAD(91) = 1.8. These parameters govern the MPI activity produced 
by JIMMY. This tune of JIMMY was arrived at by fitting the data from this analysis on the 
charged scalar particle PTsum density in the toward region for Drell-Yan production. 

In this paper the theory predictions are presented as smooth curves. These curves come from 
fits to QCD Monte-Carlo output with limited statistical accuracy. The theory curves presented 
here reproduce the QCD Monte-Carlo results (with infinite statistical accuracy) within about 2%. 

III. ANALYSIS STRATEGY 

(1) Data Sample and Event Selection 

The CDF Run II detector, in operation since 2001, is an azimuthally and forward- 
backward symmetric solenoidal particle detector [13]. It combines precision charged particle 
tracking with fast projective calorimetry and fine grained muon detection. Tracking systems are 
designed to detect charged particles and measure their momenta and displacements from the 
point of collision, termed the primary interaction vertex. The tracking system consists of a silicon 
microstrip system and an open-cell wire drift chamber, termed the Central Outer Tracker (COT) 
that surrounds the silicon. Segmented electromagnetic and hadronic sampling calorimeters 
surround the tracking system and measure the energy of interacting particles. Particles make 
showers which deposit energy and are sampled via their ionization. The muon system resides 
beyond the calorimeters. Muons are minimally ionizing particles and, hence, only deposit small 
amounts of ionization energy in the material. They are the only particles likely to penetrate both 
the tracking and five pion absorption lengths of calorimeter steel, and leave tracks in the muon 
detection system. At CDF the positive z-axis is defined to lie along the incident proton beam 
direction. The leading-jet data and lepton-pair data corresponds to an integrated luminosity of 
about 2.2 fb"' and 2.7 fb"', respectively. For both data sets we require one and only one primary 
vertex within the fiducial region |Zvertex| — 60 cm centered around the nominal CDF z —0. 

(2) Jet Selection 

Jets are selected using the MidPoint cone based algorithm with a cone size of 0.7 and 
fmerge = 0.75 [1]. For the leading-jet events we require that the highest pi jet in the calorimeter lie 
in the range |ti| < 2 or the event is rejected. 
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(3) Lepton Selection 

Dielectron events are triggered online by either one central (|ri| < 1.1) electron candidate 
with Ex > 18 GeV and a track with px > 18 GeV/c associated to it, or by two electromagnetic 
clusters with Ex > 18 GeV and |ri| < 3.2 where no track association is required. At offline level 
we consider only electrons with Et > 20 GeV and |ri| < 1 that also have a track matched to the 
calorimeter cluster. The electrons also have to pass certain quality criteria to verify that they are 
consistent with the electromagnetic shower characteristics as expected for electrons [14]. 

Dimuon events are triggered on at least one muon candidate that has a signal in one of the 
muon chambers with |ri| < land px > 18 GeV/c. The second muon candidate is not required to 
have a signal in the muon chambers but it must have hits in the COT. At offline level we 
consider only muon candidates with px > 20 GeV and |r|| < 1. All muon candidates are required 
to have calorimeter energy deposits consistent with those expected from a minimum ionizing 
particle. In addition, we employ a time-of-flight filter to remove cosmic ray muons. 

All leptons are required to be isolated from other charged fracks in the event by a distance 

of = ^(Arjf + (A^f < 0.4. 

(4) Lepton-Pair Selection 

The lepton pairs are formed by oppositely charged leptons, with the requirement that the 
z positions of the two leptons satisfy |Az | < 4 cm, to ensure that both leptons came from the same 
primary collision. For the Drell-Yan data we require that both leptons have px > 20 GeV/c and 
|ri| < 1 and that the invariant mass of the lepton-pair be in the range 70 < M(pair) < 1 10 GeV/c , 
with |ri(pair)| < 6. We chose this lepton-pair mass region because studies have shown that the 
lepton-pair backgrounds (mostly from events with QCD jets or events with a W-boson and jets) 
are negligible in the region of the Z-boson [15]. 

(5) Track Selection 

We consider charged tracks that have been measured by the central outer tracker (COT). The 
COT [16] is a cylindrical open-cell drift chamber with 96 sense wire layers grouped into eight 
altemating superlayers of stereo and axial wires. Its active volume covers 40 < r < 137 cm, 
where r is the radial coordinate in the plane transverse to the z axis, and |z| < 155 cm, thus 
providing fiducial coverage in |ri| < 1.1 to tracks originating within |z| < 60 cm. We include 
tracks in the region 0.5 < px < 150 GeV/c and |ri| < 1 where COT efficiency is high. At very 
high Px the frack resolution deteriorates. The upper limit of 150 GeV/c is chosen to prevent mis- 
measured tracks with very high px from distorting the average charged-particle density and the 
average charged-particle PTsum density. Tracks are required to hit at least two axial segments 
with more than 10 total hits and at least two stereo segments with more than 10 total hits in the 
COT. In addition, the tracks are required to point back to the primary vertex. We consider two 
track selections; loose and tight. The loose track selection requires |do| < 1.0 cm and |z - Zvtxl < 3 
cm, where do is the beam corrected transverse impact parameter and z - Zytx is the distance on the 
z-axis (beam axis) between the track and the primary vertex. The tight track selection requires 
that |do| < 0.5 cm and |z - Zvtxl < 2 cm. The loose criterion is similar to the Run 1 underlying 
event analysis [4]. 
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(6) Correcting to the Particle Level and Systematic Uncertainties 

The raw data at the detector level are corrected to the level of final state stable particles 
and are then compared with the QCD Monte-Carlo models at the generator level. The particle 
level corresponds to the event without detector efi'ects. The detector level corresponds to the 
tracks passing the good track criterion. We rely on the QCD Monte-Carlo models and the CDF 
detector simulation CDFSIM (parametrized response of the CDF II detector [17, 18]) to correct 
the measured tracks back to the stable charged particle level. Particles are considered stable if cx 
> 10 mm. The generator level charged particles have px > 0.5 GeV/c, |ri| < 1, and are kept stable 
if cx > 10 mm. Hence, to compare the corrected data with QCD Monte-Carlo model predictions 
one must keep the Kshon meson stable as well as the following baryons: A, Z, S, and Q. 

The QCD Monte-Carlo models are used to calculate the observables in Table 1.1 at the 
particle level in bins of particle jet#l pi (GEN) and at the detector level in bins of calorimeter 
jet#l Pt (CDFSIM). GEN refers to the Monte-Carlo model at the generator level and CDFSIM 
are the GEN particles after detector simulation. The detector-level data in bins of calorimeter 
jet#l Pt are corrected by multiplying by the QCD Monte-Carlo correction factor, Fcor = 
GEN/CDFSIM. This is done bin-by-bin for every observable. We refer to the ratio Fres = 
CDFSIM/GEN as the response factor for that observable with the correction factor being the 
reciprocal. Smooth curves are drawn through the QCD Monte-Carlo predictions at both the 
generator level (GEN) and the detector level (CDFSIM) to aid in comparing the theory with the 
data and also to construct the correction factors. This one step correction method simultaneously 
corrects for mis-measurement of the leading jet transverse momentum (jet energy scale) and for 
missed and/or fake tracks. 




Transverse" Charged Particle Density: SIM/GEN | 

CDF Run 2 -Leaaing Jet (pyA)- 

^atio; CDFSIM/GEN MidPoint R=0.7 |iiOet#1)|<2 



Charged Particles (|tiI<1.0, PT>0.5 GeV/c) 



50 100 150 200 250 300 350 400 450 500 
PT(jet#1 uncorrected) (GeV/c) 



"Transverse" Charged Particle Density: SIM/GEN | 



CDF Run 2 

Ratio: CDFSIM/GEN 



"Leading Jet (tiglit)" 
Midpoint R=0.7 h(jet#1)l<2 




Charged Particles {hj<1.0, PT>0.5 GeV/c) 



50 100 150 200 250 300 350 400 450 500 
PT(jet#1 uncorrected) (GeV/c) 



"Toward" Charged Particle Density: SIM/GEN | 




50 100 150 200 250 300 350 400 450 500 
PT{jet#1 uncorrected) (GeV/c) 



"Toward" Charged Particle Density: SIM/GEN | 



CDF Run 2 
Ratio: CDFSIM/GEN 



"Leading Jet (tight)" 
Midpoint R=0.7 lq(jet#1)|<2 




Charged Particles (liil<1.0, PT>0.5 GeV/c) 



50 100 150 200 250 300 350 400 450 500 
PT(jet#1 uncorrected) (GcV/c) 



Fig. 3.1. The response factors, T^.^^ = CDFSIM/GEN, for the charged-particle density, dN/drid(|), in the toward and transverse 
regions for leading-jet events. The top two plots show the response factors for PYTHIA Tune A (pyA) with tight and loose track 
cuts for the transverse (top left) and toward (top right) regions. The bottom two plots compare PYTHIA Tune A (pyA) and 
HERWIG without MPI (HW) for tight track cuts for the transverse (bottom left) and toward (bottom right) regions. The correction 
factor is the reciprocal of the response factor (T^^^ = l/F^s). 
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The correction factors are different for every observable and they are different for the 
tight and loose track selection criterion. The tight track criterion results in fewer tracks than the 
loose criterion and hence the Monte-Carlo corrections factors are different. If the Monte-Carlo 
described the data perfectly and if CDFSIM were exact, then the corrected observable would be 
identical regardless of the track selection criterion. Using PYTHIA Tune A for the leading-jet 
events and PYTHIA Tune AW for the Drell-Yan events, we find that the loose and tight track 
selections do result in nearly the same particle level result for all the observables presented in 
this analysis. The differences are used as a source of systematic error and are added in 
quadrature to the statistical errors. 

Figure 3.1 shows the response factors, Fres, for the charged-particle density in the toward 
and transverse regions for leading-jet events. The correction factors (1/Fres) are typically small 
(they differ from one by less that 10%) except in regions where the charged-particle density 
becomes large, which occurs in the toward and away regions for leading-jet production. The 
efficiency of detecting charged tracks decreases when the density of tracks becomes large. For 
this reason we restrict ourselves to the range pTGet#l) < 200 GeV/c for the toward and away 
regions, but allow the leading jet transverse momentum to extend to 400 GeV/c in the transverse 
region. For the leading-jet events we have also used HERWIG (without MP I) as well as PYTHIA 
Tune A to correct the data to the particle level. We use the differences in the corrected data as an 
additional source of systematic error (added in quadrature). For low pTGet#l) the correction 
factors become large due to the uncertainty in the jet energy scale at low energy. Also, the 
corrections from HERWIG and PYTHIA Tune A differ significantly in this region. This results in 
very large systematic errors on the data at low leading-jet transverse momentum. 

Another important effect and resultant systematic error arises from the uncertainty in the 
jet energy scale for pr of the leading jet. The CDF detector simulation does not reproduce 
perfectly the response of the calorimeters. The overall systematic uncertainty in the CDF jet 
energy scale (JES) is a function of the jet px [21]. The uncertainty is about 3% at high px and 
increases to around 8% at low px. After correcting the data to the particle level we shift PxGet#l) 
up and down by this additional uncertainty with the bin-by-bin differences in the observables in 
Table 1 . 1 used as another systematic error. The JES systematic errors are large in the toward and 
away region where the observables are varjdng rapidly with Vj(jet#l). 

We investigated the dependence of the corrected data to our upper limit of PTinax(cut) = 
150 GeV/c which was applied to all tracks. The sensitivity of the results to this choice of upper 
limit was checked by changing the upper limit to PTmax(cut) = 1 .5 x ETmax(tower). Here one 
looks, on an event-by-event basis, at all the towers in the region |ri| < 1 and sets the maximum px 
track cut to be equal to 1 .5 times the Ex of the tower with the largest transverse energy. High px 
mis-measured tracks do not deposit energy in the calorimeter. The two maximum px track cut 
methods produce slightly different correction factors, however, after correcting to the particle 
level the results are nearly identical. For the leading-jet analysis the differences were used as an 
additional systematic error. 

Although we require one and only one high quality vertex, the observables in Table 1 . 1 
can still be affected by pile-up (more than one proton-antiproton collision in the event). Tracks 
are required to point back to the primary vertex, but the track observables are affected by pile-up 
when two vertices overlap. Vertices within about 3 cm of each other merge together as one. In 
the leading-jet analysis we examined the effects of pile-up by plotting the transverse charged- 
particle density and the charged-particle PTsum density versus the instantaneous luminosity 
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(with one and only one vertex). As the instantaneous luminosity increases so does the amount of 
pile-up. We found that these observables did increase slightly with increasing luminosity 
(roughly linearly). The leading-jet observables in the transverse region are corrected for pile-up 
by extrapolating to the low luminosity limit. To correct the data, we define a low region, Linst < 
25 X 10^° cm"^s"^ (low), and a high region Linst > 25 x 10^° cm"^s"' (high), where Linst is the 
instantaneous luminosity. On a bin-by-bin basis, the ratio high/low and all/low was constructed, 
where all = high + low. The ratio high/low was found to be small (usually less than 1%) and 
could simply have been absorbed into the overall systematic errors. However, in the leading-jet 
analysis we corrected the data for pile-up by drawing a smooth curve through the ratio all/low 
and then dividing the data by this ratio. The size of the pile-up correction was then taken as the 
systematic error in making the correction and added in quadrature with the other systematic 
errors. For the Drell-Yan analysis, the pile-up corrections were less than 1% and were simply 
absorbed into the overall systematic errors. 



"Transverse" Charged Particle Density: SIIVI/GEN 



CDF Run 2 

Ratio: CDFSilVI/GEN 




"Leading jet" pyA 



a 0.8 



"Drell-Yan" pyAW 



Charged Particles (|tiI<1.0, PT>0.5 GeV/c) 



20 40 60 80 100 120 140 160 180 200 
PT(jet#1 uncorrected) or PT(lepton-pair) (GeV/c) 



Fig. 3.2. The response factors, Ty^s = CDFSIM/GEN, for the charged-particle density, dN/drid(|), in the transverse region for 
leading-jet events and for Drell-Yan events. The plots show the response factors for PYTHIA Tune A (pyA) with tight track cuts 
(leading-jet) and for PYTHIA Tune AW (pyAW) with tight track cuts (Drell-Yan). The correction factor is the reciprocal of the 
response factor (F^or = l/Fjes). 



Figure 3.2 shows the response factors, Fres = CDFSIM/GEN, for the charged-particle 
density, dN/drid(j), in the transverse region for leading-jet events and for Drell-Yan events. The 
response factors are similar, but not the same. In the Drell-Yan analysis we required the leptons 
to be isolated from other particles in the event. This biases one against a very active underlying 
event which is compensated for by the correction factor. 



IV. RESULTS 



(1) Leading-Jet and Drell-Yan Topologies 

Figure 4.1 shows the data on the density of charged particles and the scalar PTsum 
density, respectively, for the toward, away, and transverse regions for leading-jet and Drell-Yan 
events. For leading-jet events the densities are plotted as a flinction of the leading-jet px and for 
Drell-Yan events there are plotted versus the pi of the lepton pair. The data are corrected to the 
particle level and are compared with PYTHIA Tune A (leading-jet) and Tune AW (Drell-Yan) at 
the particle level. For leading-jet events at high pTOet#l) the densities in the toward and away 
regions are much larger than in the transverse region because of the toward-side and away-side 
jets. At smaU pTOet#l) the toward, away, and transverse densities become equal and go to zero 
as pTO et#l) goes to zero. If the leading jet has no transverse momentum then there can be no 
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particles anywhere. In addition, there are a lot of low transverse momentum jets and for 
PTOet#l) < 30 GeV/c the leading jet is not always the jet resulting from the hard 2-to-2 
scattering. This produces a bump in the transverse density in the range where the toward, away, 
and transverse densities become similar in size. For Drell-Yan events the toward and transverse 
densities are both small and almost equal. The away density is large due to the away-side jet. 
The toward, away, and transverse densities become equal as px of the lepton pair goes to zero, 
but unlike the leading-jet case the densities do not vanish at pT(lepton-pair) = 0. For Drell-Yan 
events with pT(lepton-pair) = the hard scale is set by the lepton-pair mass which is in the region 
of the Z-boson, whereas in leading -jet events the hard scale goes to zero as transverse 
momentum of the leading jet goes to zero. 
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Fig. 4.1. CDF data at 1.96 TeV on the density of charged particles, dN/dr|d(|), and the scalar PTsum density of charged particles, 
dPT/dr|d(|), with pj > 0.5 GeV/c and |ri| < 1 for leading-jet and Drell-Yan events as a function of the leading-jet Pt and piClepton- 
pair), respectively, for the toward, away, and transverse regions. The data are corrected to the particle level and are compared 
with PYTHIA Tune A (pyA) and Tune AW (pyAW), respectively, at the particle level. 



Figure 4.2 compares the data for leading-jet events with the data for Drell-Yan events for 
the density of charged particles and the scalar PTsum density, respectively, in the transverse 
region. The data are compared with PYTHIA Tune A (leading-jet) , Tune AW (Drell-Yan), and 
HER WIG (without MPI). For large pTGet#l) the transverse densities are similar for leading-jet 
and Drell-Yan events as one would expect. HERWIG (without MPI) does not produce enough 
activity in the transverse region for either process. HERWIG (without MPI) disagrees more with 
the transverse region of Drell-Yan events than it does with the leading-jet events. This is 
because there is no final-state radiation in Drell-Yan production so that the lack of MPI becomes 
more evident. 

Fiure 4.3 compares the data for leading-jet events with the data for Drell-Yan events for 
the average charged-particle pi and the average maximum charged-particle pi, respectively, in 
the transverse region. The data are compared with PYTHIA Tune A (leading-jet) , Tune AW 
(Drell-Yan), and HERWIG (without MPI). MPI provides a hard component to the underlying 
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event and for HER WIG (without MPI) the pi distributions in the transverse region for both 
processes are too soft, resuhing in an average pi and average PTmax that are too small. 
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Fig. 4.2. CDF data at 1.96 TeV on the density of charged particles, dN/dr|d(j), and the scalar PTsum density of charged particles, 
dPT/dr|d(|), with pj > 0.5 GeV/c and hi < 1 for leading-jet and Drell-Yan events as a function of the leading-jet py and pT(lepton- 
pair), respectively, for the transverse region. The Drell-Yan data are compared with PYTHIA Tune AW (pyAW) and the leading- 
jet data are compared with PYTHIA Tune A (pyA). Also shown are some prediction from HERWIG without MPI (HW). 
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Fig. 4.3. CDF data at 1.96 TeV on the average charged particle transverse momentum, <Pt>, and the average maximum Pt, with 
Px > 0.5 GeV/c and ir|j < 1 for leading-jet and Drell-Yan events as a function of the leading-jet pj and pT(lepton-pair), 
respectively, for the transverse region. The Drell-Yan data are compared with PYTHIA Tune AW (pyAW) and the leading-jet 
data are compared with PYTHIA Tune A (pyA). Also shown are some prediction from HER WIG without MPI (HW). 



Figure 4.4 compares the data for leading-jet events with the data for Drell-Yan events for 
the density of charged particles and the scalar PTsum density, respectively, for the transMAX 
and transMIN regions. The data are compared with PYTHIA Tune A (leading-jet) , Tune AW 
(Drell-Yan), and HERWIG (without MPI). For events with large initial-state or final-state 
radiation the transMAX region would contain the third jet in high-px jet production or the second 
jet in Drell-Yan production. Thus, the transMIN region is very sensitive to the modeling of the 
underlying event. 
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Fig. 4.4. CDF data at 1.96 TeV on the density of charged particles, dN/dr|d(|), and the scalar PTsum density of charged particles, 
dPT/dr|d(|), with p-r > 0.5 GeV/c and jril < 1 for leading-jet and Drell-Yan events as a fianction of the leading-jet py and pT(lepton- 
pair), respectively, for the transMAX and transMESf regions. The Drell-Yan data are compared with PYTHIA Tune AW (pyAW) 
and the leading-jet data are compared with PYTHIA Tune A (pyA). Also shown are some prediction from HERWIG without MPI 
(HW). 



Figure 4.5 compares the data for leading-jet events with the data for Drell-Yan events for 
the density of charged particles and the scalar PTsum density for transDIF = transMAX - 
transMIN. The data are compared with PYTHIA Tune A (leading-jet) and Tune AW (Drell-Yan). 
The transDIF region is sensitive to the hard initial and final-state radiation and is predicted to be 
very similar in the two processes. Fig. 4.5 also compares the data for leading-jet events with the 
data for Drell-Yan events for the density of charged particles and the scalar PTsum density in the 
away region. The away-side jet pseudorapidity distribution and type (quark or gluon) is different 
for leading-jet and Drell-Yan events so we do not expect the away region to be the same and it is 
not. However, PYTHIA Tune A and Tune AW describe the data very well. 
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Fig. 4.5. CDF data at 1.96 TeV on the density of charged particles, dN/dr|d(|), and the scalar PTsum density of charged particles, 
dPT/dr|d(|), with pj > 0.5 GeV/c and |ri| < 1 for leading-jet and Drell-Yan events as a fianction of the leading-jet Pt and pjClepton- 
pair), respectively, for the transDIF region (transDIF = transMAX - transMESf) and the away region. The Drell-Yan data are 
compared with PYTHIA Tune AW (pyAW) and the leading-jet data are compared with PYTHIA Tune A (pyA). 
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Fig. 4.6. CDF data at 1.96 TeV on the density of charged particles, dN/drid(|), the scalar PTsum density of charged particles, 
dPT/drid(|), and the average charged-particle px, with pj > 0.5 GeV/c and |ri| < 1 for Drell-Yan events as a function of pT(lepton- 
pair) for the toward, transverse, and transMIN regions compared with PYTHIA Tune AW (pyAW). 
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(2) The Underlying Event in Drell-Yan Production 



Figure 4.6 compares the data in the toward region with the data in the transverse region 
for Drell-Yan events for the density of charged particles, the scalar PTsum density, and the 
average charged-particle px. The data are compared with PYTHIA Tune AW. For high transverse 
momentum lepton-pair production, particles from initial-state radiation are more likely to 
populate the transverse region than the toward region and hence the densities are slightly larger 
in the transverse region. PYTHIA Tune AW describes this very nicely. 
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Fig. 4.7. CDF data at 1.96 TeV on the density of charged particles, dN/drid(|), the scalar PTsum density of charged particles, 
dPT/dr|d(|), and the average charged-particle pj, and the average maximum charged-particle pj. with pj > 0.5 GeV/c and jr|j < 1 
for Drell-Yan events as a function of pT(lepton-pair) for the toward and the transMIN regions. The data are compared with 
HERWIG without MPI (HW), HERWIG with JIMMY MPI (JIM), and three PYTHIA Tunes (pyAW, pyDW, ATLAS). 



The most sensitive regions to the underlying event in Drell-Yan production are the 
toward and the transMIN regions, since these regions are less likely to receive contributions from 
the away-side jet and from initial-state radiation. Fig. 4.7 shows the data for Drell-Yan events 
for the density of charged particles and the scalar PTsum density, respectively, in the toward and 
transMIN regions. The data are compared with PYTHIA Tune AW, Tune DW, the PYTHIA 
ATLAS tune, HERWIG (without MPI), and HERWIG (with JIMMY MPI). The densities are smaller 
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in the transMIN region than in the toward region and this is described weU by PYTHIA Tune AW. 
Comparing HERWIG (without MPI) with HERWIG (with JIMMY MPI) clearly shows the 
importance of MPI in these regions. Tune AW and Tune DW are very similar. The ATLAS 
tune and HERWIG (with JIMMY MPI) agree with Tune AW for the scalar PTsum density in the 
toward and transMIN regions. However, both the ATLAS tune and HERWIG (with JIMMY MPI) 
produce too much charged-particle density in these regions. The ATLAS tune and HERWIG 
(with JIMMY MPI) fit the PTsum density, but they do so by producing too many charged 
particles. They both have too soft a px spectrum in these regions. This can be seen clearly in Fig. 
4.7 which shows the data for Drell-Yan events on the average charged-particle px and the 
average maximum charged-particle px, in the toward region compared with the QCD Monte- 
Carlo models. 
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Fig. 4.8. CDF data on the density of charged particles, dN/drid(|), with p-p > 0.5 GeV/c and [r|[ < 1 for Drell-Yan events in the 
toward region and for leading-jet events in the transverse region at 1 .96 TeV (Tevatron) together with the predictions of PYTHIA 
Tune DW (pyDW) and Tune DWT (pyDWT) at the 14 TeV (LHC). Tune DW and Tune DWT are identical at 1.96 TeV 
(Tevatron). Also shown are CDF data at 1.96 TeV (Tevatron) on the average charged particle transverse momentum, <Pt>, for 
Drell-Yan events in the toward region compared with HERWIG without MPI (HW) and PYTHIA Tune DW (pyDW) together with 
extrapolations to the 14 TeV (LHC). 



(3) Extrapolating to the LHC 

Figure 4.8 shows the extrapolation of PYTHIA Tune DW, Tune DWT, and HERWIG (without 
MPI) to 14 TeV (LHC) for the density of charged particles and the average transverse 
momentum of charged particles with px > 0.5 GeV/c and |ri| < 1 in the towards region of Drell- 
Yan production. The underlying event activity is the same for proton-proton and proton- 
antiproton collisions. For HERWIG (without MPI) the toward region of Drell-Yan production 
does not change much in going from the Tevatron to the LHC. Fig. 4.8 also shows the 
extrapolation of PYTHIA Tune DW and Tune DWT to 14 TeV (LHC) for the transverse density 
of charged particles (pt > 0.5 GeV/c, |ri| < 1) for leading-jet events. Models with multiple- 
parton interactions like PYTHIA Tune DW and Tune DWT predict that the underlying event will 
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become much more active (with larger <Pt>) at the LHC. PYTHIA Tune DW predicts about a 
factor of two increase in the activity of the underlying event as measured by the charged-particle 
density (pj > 0.5 GeV/c, |ri| < 1) in the towards region of Drell-Yan production and the 
transverse region in leading-jet events. Tune DWT used the default value for PARP(90) and 
predicts an even greater increase in the activity of the underlying event at the LHC. However, 
Tune DWT produces less activity than Tune DW in the underlying event at energies below 1.96 
TeV and the CDF data at 630 GeV [9] favor Tune DW over Tune DWT. 

(4) <Pt> versus the Multiplicity: Min-Bias and Drell-Yan Events 

The total proton-antiproton cross section is the sum of the elastic and inelastic components, 
CTtot = ctel + CTiN. The inelastic cross section consists of three terms; single diffraction, double- 
diffraction, and everything else (referred to as the hard core), aiN = osd + cjdd + cjhc- For elastic 
scattering neither of the beam particles breaks apart. For single and double diffraction one or 
both of the beam particles are excited into a high mass color singlet state (i.e. N states) which 
then decays. Single and double diffraction also corresponds to color singlet exchange between 
the beam hadrons. When color is exchanged, the outgoing remnants are no longer color singlets 
and one has a separation of color resulting in a multitude of quark-antiquark pairs being pulled 
out of the vacuum. The hard core component, ohc, involves color exchange and the separation 
of color. However, the hard core contribution has both a soft and hard component. Most of the 
time the color exchange between partons in the beam hadrons occurs through a soft interaction 
with no high transverse momentum and the two beam hadrons ooze through each other 
producing lots of soft particles with a uniform distribution in rapidity and many particles flying 
down the beam pipe. Occasionally there is a hard scattering among the constituent partons 
producing outgoing particles and jets with high transverse momentum. 

Minimum bias (min-bias) is a generic term which refers to events that are selected with a 
loose trigger that accepts a large fraction of the inelastic cross section. All triggers produce some 
bias and the term min-bias is meaningless until one specifies the precise trigger used to collect 
the data. The CDF min-bias trigger consists of requiring at least one charged particle in the 
forward region 3.2 < t] < 5.9 and simultaneously at least one charged particle in the backward 
region -5.9 <r\< -3.2. Monte-Carlo studies show that the CDF min-bias trigger collects most of 
the Ohc contribution plus small amounts of single and double diffraction [20]. 

Minimum bias collisions are a mixture of hard processes (perturbative QCD) and soft 
processes (non-perturbative QCD) and are, hence, very difficult to simulate. Min-bias collisions 
contain soft beam-beam remnants, hard QCD 2-to-2 parton-parton scattering, and multiple 
parton interactions (soft & hard). To correctly simulate min-bias collisions one must have the 
correct mixture of hard and soft processes together with a good model of the multiple-parton 
interactions. We have seen that multiple parton interactions are a significant component of the 
underlying event in high px jet production and in Drell-Yan lep ton-pair production. Multiple- 
parton interactions are also an important component in min-bias collisions. Min-bias collisions 
are not the same as the underlying event in a hard-scattering process, since the rate at which MPI 
occurs is different, but they are related. In selecting a hard-scattering process such as high px jet 
production or in selecting lepton-pair in the mass region of the Z-boson corresponds to selecting 
a small fraction of min-bias collisions that are very central; the initial proton and antiproton 
collide with small impact parameter. For these central collisions the probability of additional 
parton-parton collisions is higher than it is for an average min-bias event. 
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Fig. 4.9. (top) CDF Min-Bias data on the average pj of charged particles versus the multiplicity for charged particles with pj > 
0.4 GeV/c and Inl < 1 from Ref [20]. The data are compared with PYTHIA Tune A (pyA), the PYTHIA ATLAS tune, and 
PYTHIA Tune A without MPI (pyAnoMPI). {middle) CDF data on the average px of charged particles versus the multiplicity for 
charged particles with px > 0.5 GeV/c and jr|| < 1 for Drell-Yan events, (bottom) CDF data on the average Pt of the lepton-pair 
versus the multiplicity for charged particles with pj > 0.5 GeV/c and [r|[ < 1. The Drell-Yan data are compared with PYTHIA 
Tune AW, the PYTHIA ATLAS tune, HERWIG without MPI (HW), and HERWIG with JIMMY MPI (JIM). 



The first model that roughly described min-bias collisions at CDF was PYTHIA Tune A. 
However, Tune A was not tuned to fit min-bias collisions. It was tuned to fit the activity in the 
underlying event in high transverse momentum jet production [4]. However, PYTHIA uses the 
same px cut-off for the primary hard 2-to-2 parton-parton scattering and for additional multiple 
parton interactions (MPI). Hence, fixing the amount of multiple parton interactions by setting 
the Pt cut-off allows one to run the hard 2-to-2 parton-parton scattering all the way down to 
PT(hard) = without hitting a divergence. For PYTHIA the amount of hard scattering in min-bias 
is, therefore, related to the activity of the underlying event in hard-scattering processes. Neither 
HERWIG (without MPI) or HERWIG (with JIMMY MPI) can be used to describe min-bias events 
since they diverge as pT(hard) goes to zero. 
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Figure 4.9 shows CDF min-bias data corrected to the particle level at 1.96 TeV on the average 
Pt of charged particles, <Pt>, versus the multiplicity for charged particles with px > 0.4 GeV/c 
and |ri| < 1 from Ref [20]. The data are compared with PYTHIA Tune A, the PYTHIA ATLAS 
tune, and PYTHIA Tune A without MPI (pyAnoMPI). The average px is an important observable. 
The rate of change of <px> versus charged multiplicity is a measure of the amount of hard versus 
soft processes contributing to min-bias collisions and it is sensitive the modeling of the multiple- 
parton interactions [21]. If only the soft beam-beam remnants contributed to min-bias collisions 
then <px> would not depend on charged multiplicity. If one has two processes contributing, one 
soft (beam-beam remnants) and one hard (hard 2-to-2 parton-parton scattering), then demanding 
large multiplicity will preferentially select the hard process and lead to a high <Pt>. However, 
we see that with only these two processes <px> increases much too rapidly as a function of 
multiplicity (see pyAnoMPI). Multiple-parton interactions provide another mechanism for 
producing large multiplicities that are harder than the beam-beam remnants, but not as hard as 
the primary 2-to-2 hard scattering. PYTHIA Tune A gives a fairly good description of the <Pt> 
versus multiplicity, although not perfect. PYTHIA Tune A does a better job describing the data 
than the ATLAS tune. Both Tune A and the ATLAS tune include multiple-parton interactions, 
but with different choices for the color connections [22]. 
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Fig. 4.10. (top) CDF data at 1.96 TeV on the average Pt of charged particles versus the multiplicity for charged particles with pj 
> 0.5 GeV/c and [r|[ < 1 for Drell-Yan events in which px(lepton-pair) < 10 GeV/c. The data are compared with PYTHIA Tune 
AW (pyAW), the PYTHIA ATLAS tune, HERWIG without MPI (HW), and HERWIG with JIMMY MPI (JIM), {bottom) 
Comparison of the average px of charged particles versus the charged multiplicity for Min-Bias events from Ref. 20 with the 
Drell-Yan events with pT(lepton-pair) < 10 GeV/c from this analysis. The Min-Bias data require pj > 0.4 GeV/c and are 
compared with PYTHIA Tune A (pyA), while the Drell-Yan data require pj > 0.5 GeV/c and are compared with PYTHIA Tune 
AW (pyAW). 

Figure 4.9 also shows the data at 1.96 TeV on the average px of charged particles versus the 
multiplicity for charged particles with px > 0.5 GeV/c and |ti| < 1 for Drell-Yan events from this 
analysis. HERWIG (without MPI) predicts the <px> to rise too rapidly as the multiplicity 
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increases. This is similar to the pyAnoMPI behavior in min-bias collisions. For HERWIG 
(without MPI) large multiplicities come from events with a high px lepton-pair and hence a large 
Pt away-side jet. This can be seen clearly in Fig. 4.9 which also shows the average pi of the 
lepton-pair versus the charged multiplicity. Without MPI the only way of getting large 
multiplicity is with high-pT(lepton-pair) events. For the models with MPI one can get large 
multiplicity either from high-pxClepton-pair) events or from MPI and hence <PT(lepton-pair)> 
does not rise as sharply with multiplicity in accord with the data. PYTHIA Tune AW describes 
the Drell-Yan data fairly well. 

Figure 4.10 shows the data at 1.96 TeV on the average px of charged particles versus the 
multiplicity for charged particles with px > 0.5 GeV/c and Itll < 1 for Drell-Yan events in which 
Px(lepton-pair) < 10 GeV/c. We see that <px> still increases as the multiplicity increases 
although not as fast. If we require px(lep ton-pair) < 10 GeV/c, then HERWIG (without MPI) 
predicts that the <px> decreases slightly as the multiplicity increases. This is because without 
MPI and without the high px away-side jet which is suppressed by requiring low px of the lepton 
pair, large multiplicities come from events with a lot of initial-state radiation and the particles 
coming from initial-state radiation are soft. PYTHIA Tune AW describes the behavior of <px> 
versus the multiplicity fairly well even when we select px(lepton-pair) < 10 GeV/c. 

Figure 4.10 also shows a comparison of the average px of charged particles versus the charged 
multiplicity for min-bias events from Ref 20 with the Drell-Yan events with px(lepton-pair) < 10 
GeV/c. There is no reason for the min-bias data to agree with the Drell-Yan events with 
Px(lepton-pair) < 10 GeV/c. However, they are remarkably similar and described fairly well by 
PYTHIA Tune A and Tune AW, respectively. This strongly suggests that MPI are playing an 
important role in both these processes. 

V. Summary & Conclusions 

Observables that are sensitive to the underlying event in high transverse momentum jet 
production (leading-jet events) and Drell-Yan lepton pair production in the mass region of the Z- 
boson (Drell-Yan events) have been presented and compared with several QCD Monte-Carlo 
model tunes. The data are corrected to the particle level and compared with the Monte-Carlo 
models at the particle level. The underlying event is similar for leading-jet and Drell-Yan events 
as one would expect. This analysis provides data that can be used to test and improve the QCD 
Monte-Carlo models of the underlying event that are used to simulate hadron-hadron collisions. 
The data presented here are also important for tuning the new QCD Monte-Carlo multiple-parton 
interaction (MPI) models [21, 20]. 

PYTHIA Tune A and Tune AW do a good job in describing the data on the underlying 
event observables for leading-jet and Drell-Yan events, respectively, although the agreement 
between theory and data is not perfect. The leading-jet data show slightly more activity in the 
underlying event than PYTHIA Tune A. PYTHIA Tune AW is essentially identical to Tune A for 
leading-jet events. All the tunes MPI agree better than HERWIG without MPI. This is especially 
true in the toward region in Drell-Yan production. Adding JIMMY MPI to HERWIG greatly 
improves the agreement with data, but HERWIG with JIMMY MPI produces a charged-particle pi 
spectrum that is considerably softer than the data. The PYTHIA ATLAS tune also produces a 
charged-particle px spectra that is considerably softer than the data. 

The behavior of the average charged-particle px versus the charged-particle multiplicity is 
important. The rate of change of <px> versus charged multiplicity is a measure of the amount of 
hard versus soft processes contributing, and it is sensitive the modeling of the multiple-parton 
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interactions. PYTHIA Tune A and Tune AW do a good job in describing the data on <pt> versus 
multiplicity for min-bias and Drell-Yan events, respectively, although again the agreement 
between theory and data is not perfect. The behavior of <Pt> versus multiplicity is remarkably 
similar for min-bias events and Drell-Yan events with pT(Iepton-pair) < 10 GeV/c, suggesting 
that MPI are playing an important role in both these processes. 

Models with multiple-parton interactions like PYTHIA Tune DW predict that the 
underlying event will become much more active (with larger <Pt>) at the LHC. For HERWIG 
(without MPI) the toward region of Drell-Yan production does not change much in going from 
the Tevatron to the LHC. PYTHIA Tune DW predicts about a factor of two increase in the 
activity of the underlying event in going from the Tevatron to the LHC as measured by the 
charged-particle density (px > 0.5 GeV/c, |ri| < 1) in the towards region of Drell-Yan production 
and the transverse region in leading-jet events. Tune DWT predicts an even greater increase in 
the activity of the underlying event at the LHC. However, Tune DWT produces less activity 
than Tune DW in the underlying event at energies below 1 .96 TeV. Tune DW does a better job 
in fitting the CDF underlying event data at 630 GeV [9], and is hence favored over Tune DWT. 
At present, PYTHIA tunes with PARP(90) around the value of Tune AW and Tune DW (~ 0.25) 
seem to be preferred. We will learn a lot about the energy dependence of MPI by comparing the 
Tevatron results with the early LHC measurements and precise measurements at the LHC require 
good modeling of the underljdng event. 
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